arXiv: 1505.06673vl [astro-ph.HE] 25 May 2015 


Astronomy <fe Astrophysics manuscript no. OscillatingTorusVS @ ESO 2015 

May 26, 2015 


Twin peak HF QPOs as a spectral imprint of dual oscillation 

modes of accretion tori 

P. Bakala^, K. Goluchova^, G. Torok^, E. Sramkova^, M. A. Abramowicz^’^’^, F. H. Vincent^, and G. P. 

Mazur^’^ 

^ Institute of Physics, Faculty of Philosophy and Science, Silesian University in Opava, Bezrucovo nam. 13, CZ-746 01 
Opava, Czech Republic 
e-mail: pavel.bakala@fpf.slu.cz 
e-mail: katka.g@seznam.cz 
e-mail: gabriel.torok@gmail.com 

^ Nicolaus Copernicus Astronomical Center, ul. Bartycka 18, PL-00-716 Warszawa, Poland 
e-mail: fvincent@camk.edu.pl 

® Institute of Physics, Polish Academy of Sciences, aleja Lotnikow 32/46, PL-02-668 Warszawa, Poland 
e-mail: gmazur@camk. edu. pi 

Physics Department, Gothenburg University, SE-412-96 Goteborg, Sweden 
e-mail: marek. abramowicz@physics .gu.se 

Received ; accepted 


Abstract 


Context. High frequency (millisecond) quasi-periodic oscillations (HF QPOs) are observed in the X-ray power-density 
spectra of several microquasars and low mass X-ray binaries. Two distinct QPO peaks, so-called twin peak QPOs, 
are often detected simultaneously exhibiting their frequency ratio close or equal to 3/2. A widely discussed class of 
proposed QPOs models is based on oscillations of accretion toroidal structures orbiting in the close vicinity of black 
holes or neutron stars. 

Aims. Following the analytic theory and previous studies of observable spectral signatures, we aim to model the twin 
peak QPOs as a spectral imprint of specific dual oscillation regime dehned by a combination of the lowest radial 
and vertical oscillation mode of slender tori. We consider the model of an optically thick slender accretion torus with 
constant specific angular momentum. We examined power spectra and fluorescent Ka iron line profiles for two different 
simulation setups with the mode frequency relations corresponding to the epicyclic resonance HF QPOs model and 
modified relativistic precession QPOs model. 

Methods. We use relativistic ray-tracing implemented in parallel simulation code LSDplus. In the background of the 
Kerr spacetime geometry, we analyze the influence of the distant observer inclination and the spin of the central compact 
object. Relativistic optical projection of the oscillating slender torus is illustrated by images in false colours related to 
the frequency shift. 

Results. We show that performed simulations yield power spectra with the pair of dominant peaks corresponding to the 
frequencies of radial and vertical oscillation modes with the proper ratio equal to 3/2 on a wide range of inclinations 
and spin values. We also discuss exceptional cases of a very small and very high inclination as well as unstable high spin 
relativistic precession-like configuration predicting constant frequency ratio equal to 1/2. We demonstrate significant 
dependency of broadened Ka iron line profiles on the inclination of the distant observer. 

Conclusions. This study presents a further step towards the proper model of oscillating accretion tori producing HF 
QPOs. More realistic future simulations should be based on incorporation of the resonant coupling of oscillation modes, 
influence of torus opacity and pressure effects on the mode frequencies and the torus shape. 

Key words. Accretion, accretion disks - Black hole physics - Relativistic processes 


1. Introduction 

High frequency quasi-periodic oscillations (HF QPOs) have 
been observed in several microquasars and low mass X- 
ray binaries (LMXBs). Their frequencies are roughly com¬ 
parable to frequencies of orbital motion of test particles 
in the vicinity of the central compact object. The black 
hole HF QPOs occur at frequencies that are characteristic 
for a particular source and seem to be constant in time. 
Strictly speaking, the current observational data enable to 
distinguish HF QPO s only for fairly specific spectral states 
(|Belloni et al.l[201^ . However, provided that HF QPO with 
lower amplitudes and different frequencies exist in the re¬ 


maining time, their detection goes beyond our present tech¬ 
nological capabilities. 

Two kinds of sharp HF QPO peaks ca n be distinguished , 
so-called lower and upper HF Q POs (see Ivan der Klisl[2004 
iRemillard &: McClintockI l2006l for a review). If both HF 
QPOs peaks are observed simultaneously (twin-peak HF 
QPOs), the ratio of frequencies of upper and lower HF QPO 
peaks is often close to 3:2 indicati ng the possible presence of 
unspecihed resonant ph enomena (jAbramowicz &: Kluznia^ 
l2001l : lT6rok et al.]l2005D n. At present, there is no consensus 

^ In the case of LMXBs, there are indications that the ratio of 
twin peak HF QPOs frequencies is clustered not only around 3:2, 
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on the QPO nature. However, the inverse mass scali ng of 
the QPOs frequencies (|Abraniowicz fc Kluznialjl2001[l pro¬ 
vides a strong argument for interpreting the observed QPO 
peaks by the frequencies of perturbed orbital motion in the 
strong gravitational field or oscillation of some accretion 
structures, accretion structures. Such interpreting natu¬ 
rally promises an attractive possi bility to measure the mass 


20011 lAbramowicz & Frae.ildl201^ Ineram & Mottall20l4 

Motta et al.ll20l4 iTbrokIboOfil: iTbrok et al.l 2012fl. Manv 

QPOs models were proposed (e.g. lAloar fc ShahamI 

1985; 

Lamb et al. 19851 Miller et al. 19981: Psaltis et al. 

1999; 

Stella fc Vietril 19991: Abramowicz fc Kluzniaki 2001. 

2004; 

Katol 2 OOII: luitarchuk fc WoodI 12002: IRezzolla et al. 

2003; 

Schnittman & Bertschinger 2004: lPetrill2005l: Zhane 

2005; 

Tagger fc Varnierd 20061: Katol 2007: Stuchlfk et al. 

2008; 

MukhoDadhvavl 20091 

Cadez et al.1 20081: iKostic et al. 

2009; 

Germana et al. 

2009 

: Germanal 120131: iLai et al. 

2013; 

Pechacek et al. 

120081 120131 and others), however 

each 


model still faces several difficulties. Moreover, the capa¬ 
bilities of the present X-ray observatories (e.g., Rossi X-ray 
Timing Explorer - RTXE) are insufficient to adequately an¬ 
alyze the harmonic content of the power spectra of observed 
lightcurves, which can be crucial to distinguishing between 
parti cular QPO models (jBakala et al.l I20l4 iKaras et al.l 
I 2 OI 4 II . Hopefully, the proposed future instruments repre¬ 
sented b y e.g. Large Obser vatory for X-ray Timing (LOFT) 
project (IFeroci et al.ll2014ll . targeted to explore strong grav¬ 
ity environment, will advance QPOs observational possibil¬ 
ities. 


A specific class of QPO models assumes oscillations ex¬ 
cited in accretion tori. The first QPO model involving nu- 
merically mod e lled t hick accretion tori was developed by 
iRezzolla et al.l (|2003ll and related light curves and powe r 
spectra were analyzed by ISchnittman fc Rezzollal (l2006ll . 
Light curves and power spectra of radially and vertically os¬ 
cillating slender torus with a circular cross-section were in- 
ves tigated in th e backg round of the Schwarzschild geometry 
by iBursa et al.l (12004 1. Numerical simulations of epicyclic 
modes of t ori oscillat i ons w ere compared to the analytical 
results by ISramkov^ (|2005ll . Later on, there has arose of 
studies devoted to mor e realistic analytic treatment of oscil¬ 
lating slender tori le.g. lAbramowicz et al.ll200'^lBlaes et al.l 
1200611 . Such studies were extended to the case of non-slender 
tori bv IStraub fc Sramkoval (|2009ll . 

In this paper, we examine timing properties of the 
numerically simulated flux emitted from a slender, poly¬ 
tropic, perfect fluid, non-self-gravitating accretion torus 
with constant specific angular momentum, which oscillates 
simultaneously in radial and vertical directions. 0 The ar ¬ 
ticle is a follow-up of th e studies of iMazur et al.l (|2013ll 
and IVincent et al.l (|2014ll who investigated the observable 
signatures of simple time-periodic slender torus deforma¬ 
tions as well as slender torus m otion described by the set 
of oscillation modes derived bv iBlaes et al.l (|2006ll . Using 
their results, we took the next step towards a fully re¬ 


but also 4:3 and 5:4 ratios are observed (see, e.g.. iTorok et all 
120071.12008allbira') . 

^ The slender torus geometry considered here is in good ac - 
cord with the truncated disk m odel (see, e.g., iDone et ^1.1120071 1. 
Also, the power spectral fits of lingram fc Pond (|2012al i predict 
that the hot inner flow (corona) has a small scale height in the 
relevant spectral state. 


alistic model of HF QPO based on oscillations of accre¬ 
tion tori. In order to model simultaneously detected twin 
peaks HF QPOs pairs, we define a new dual oscillation 
mode as a linear combination of two lowest slender torus 
oscillation modes: radial and vertical ones. We assume that 
the observed twin peaks HF QPOs can be identified with 
the pairs of the most prominent peaks in the modelled 
power spectra. We examine two different setups of the 
dual oscillation regime corresponding to the twin peaks HF 
QPOs frequency relation of epicyclic r esonance HF QPOs 
model (jAbramowicz fc Kluzniaki l200ltl and slightly mod¬ 
ified a nalogous relations of relativistic prece ssion QPOs 
model (jStella fc Vietrilfr999l: iTdrbk et al.ll2012[ ) . Those two 
competing QPO s models are prob ably the most discussed 
ones at present (jFeroci et al.ll20l3l . 

Our simulations are performed in the background of the 
Kerr geometry corresponding to the case of microquasars 
with the black hole binary component. It was shown that 
the spacetime around slowly rotating high-mass neutron 
stars can be fairly well approximated by the Kerr metric 
(|T5rok et al.l 1201(1120121 : lUrbanec et al.f[20T^ . Therefore, 
such a finding sets conditions and constrains for the appli¬ 
cability of presented results for the case of LMXBs with 
neutron star. 

The commonly accepted model of X-ray energy spec¬ 
trum of microquasars and LMXBs assumes that the illumi¬ 
nation of the cold accretion disk or torus by the primary 
component of X-ray spectrum, interpreted as the inverse 
Compton scattering of thermal photons in a hot corona, 
produces spectral lines by fluorescence. The strongest ob¬ 
served line is the Ka iron line located at 6.4 keV in the rest 
frame. Observed broad profiles of the spectral lines are sub¬ 
stantially influenced by the spacetime metric, the geometry 
of the emitting region and the distant observer incl i natio n 
([Fabian et al.l Il98ft ICadez fc Calvanil 120051 : iBambil 120131) . 
In order to develop an additional tool for distinguishing 
various configurations of radiating slender torus, we com¬ 
pute related Ka iron line profiles. 

The article consits of the following parts: Section[5]is de¬ 
voted to the theory which describes the slender torus model 
and its dual oscillation regime. Section [3] provides more de¬ 
tails of the investigated particular setup of the model of 
slender torus. Section [4] describes our numerical implemen¬ 
tation of relativistic ray-tracing and the following construc¬ 
tion of lightcurves, power spectra and iron Ka line profiles. 
Section [5] is devoted to the methodology of simulations per¬ 
formance and results analysis. Section [5] shows the obtained 
results, and Section [7] gives conclusions and future research 
perspectives 


2. Slender torus 

2.1. Equilibrium torus configuration 

We consider an axisymmetric, non self-gravitating, perfect 
fluid, constant specific angular momentum, circularly or¬ 
biting accretion torus in the background of the Kerr geom¬ 
etry. Using the (— I- -|— 1-) signature and geometrical units 
(c = G = M = 1), the line element of the Kerr spacetime 
in Boyer-Lindquist coordinates parameterized by specific 


2 















































































































































Bakala et al.: Twin peak HF QPO as a imprint of dual oscillation modes 


angular momentum (spin) a reads 


, .. , 2 r 

ds^ = - 1 - — 


Edr 


9 1 ^i9 

dt^ -— sm^ 9dtd(p + -^ dr^ 

l \ 

y 2ra^ sin^ 9'\ , , ^ , 9 
a H-- I sm'* 9 dip , 


( 1 ) 


where E = cos^ 9 and A = — 2r + a^. Moreover, 

we assume that the radial extent of the torus cross-section 
is small compared to its central radius. The perfect fluid 
forming the torus is described by polytropic equation of 
state (with a polytropic constant K and polytropic index 
n). In such a case, the energy density e is a function of 
the pressure p and the mass density p. The pressure is a 
function of mass density only. The torus surface is given 
by the zero pressure, while the pressure gradient is equal 
to zero in th e torus centre. Using the conservation law, one 
can get Isee lAbramowicz et al.ll200^ iBlaes et al.ll200ffl 


P Po 

- = — fir,9), 
P Po 


X — (-yf 




Pro 


Pro 


P" = 


2nclo 


( 4 ) 


where Mq denotes the four-velocity of the fluid and square 
of the sound speed c^q << 1 can be expressed in the form 

_ n -I- I po 

.dpJo 


= 


( 5 ) 


fix,y) = l-uj^x -Ugy . 


where 



Log = J 1 - 


( 9 ) 


( 10 ) 


2.2. Oscillating slender torus 

Let us assume small pressure perturbation in the form 

^poc (11) 

where m is azimuthal wave number (m € N) and fl is os¬ 
cillation angular frequency. It is useful to introduce a new 
variable - eigenfuction Wi , which is the function of pertur¬ 
bation in pressure in the form 


W, = - 


Sp 


where related eigenfrequency ai reads 

(7^ = 


( 12 ) 


(13) 


The perturbative surface equation /( r, 9) and torus four- 
velocity can be expressed in the form (| Vincent et ahll^ldfl 


(2) fir,9) = fir,9)- 


1 Po„,t 


where the surface function / is constant at isobaric and iso¬ 
density surfaces. Naturally, the surface function vanishes at 
the torus surface. Here and below, values of the quantities 
denoted with a subscript q are taken in the centre of the 
torus in the equilibrium state. 

Following lAbramowicz et id] (j2006ll , we use the new ra¬ 
dial and vertical coordinates 


= Un + Re . . 

\ Po + eo V Pxr- 


ulRe{Wi}ai, 
ipo (dWi 


n+\po ° 


)}■ 


(14) 

(15) 


with zero in the torus centre. Parameter P determines the 
torus thickness and is given by formula 


n po 

Keplerian angular velocity OLko in the Kerr spacetime reads 
O.K 0 = l^/{ry'^ + a). ( 6 ) 

In the case of slender torus, the parameter P must fulfil the 
condition 

/3 « 1 ■ ( 7 ) 

In such a coordinate frame, the surface function /(r, 9) can 
be rewritten as 


Discrete set of eigefunctions Wi and related eigenfrequen- 
cies Oj desc r ibing different oscillation modes was found by 
iBlaes et al.l (|2006ll solving Papaloizou-Pringle equation for 
slender torus case corresponding to the condition 0 —>• 0 
(jAbramowicz et a02006tlPapaloizou fc Prin^ll984ll . Here 
we use only two simplest solutions: radial oscillation mode 
and vertical oscillation mode. In the case of radial oscilla¬ 
tion mode, the eigenfuction and related eigenfrequency are 
given as 

Wr = , Or = cuMko , (16) 

where Or is a free parameter related to the amplitude of 
oscillations. Then the surface equation m takes the form 

1 — — u)gy^ — ArXCOs{mrP — D^t) = 0, (17) 

where amplitude Ar is given by the formula 


A — ^ Po t 

-/i-T* — ^ Uq(T'p(X'p ^ 


(18) 


n + lpo 

and mode oscillation angular frequency reads 

= Or + flKorrir . (19) 

In order to rewrite the surface equation (HU to the form 
1 — uj'lix + Sx)^ — Ogy^ = 0 , ( 20 ) 


(8) the term 


-f- cos(mrP — firt) 


( 21 ) 


must be added. Then the displacement Sx can be expressed 
by the formula 


Sx = cos {rUrP — flrt) . 
2u}~ 


( 22 ) 


are radial and vertical epicyclic frequencies of free test par- 
ticles in the centre of the torus scaled to Djfn (see, e.g., 

Aliev &: GaltsovIlOSltlAbramowicz &: Kluzniakl200^lAliev| 

2008|1 . We can easily see, that the slender torus cross-section 
(given by /(x, y) = 0) has elliptical shape in the x-y plane. 


In this approximation, the added term should be small 
enough, which corresponds to the condition for radial am¬ 
plitudes 

§ « 1- P3) 
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The radial component of the surface four-velocity is simply 
given by a derivative of the displacement 5x with respect 
to proper time r as 

r- _ ^ _ dx /3ro _ d{-5x) /grp 
~ dr ~ dr ~ dr ’ 

and it takes the covariant form 

A 

uSur ^ (v^)o TT^^KOul Sm^mr^p - ^rt) ■ (25) 

Analogously, in the case of the vertical oscillation mode, 
the eigenfunction and related eigenfrequency are given as 

Wg = agyed'^^^-^o^^ , ag = UJg^KO ■ (26) 

Using the condition for vertical amplitudes << 1, the 
surface equation, displacement 5y and surface four-velocity 
component u™’' for the vertical mode can be expressed as 

l-ulx"^ -ul[y + 5yf = Q, (27) 

dy = cos {mgip - Qgt) , (28) 

ZtOg 

A-ff , 

= /3ro -;—^KoUo shi{mgip - ^gt ), (29) 

ZiJg 

where amplitude Ag given by formula 

Ag = — A—ulagag , (30) 

n -I- 1 po 

and mode oscillation angular frequency reads 

^}g = ag + ^Komg ■ (31) 

2.3. Dual oscillation mode 

In such an approximation, radial and vertical oscillations 
are independent, and we can easily combine them into a 
new surface equation 

1-ujI{x + dxf' -ujl{y + Syf = 0. (32) 

The equation above defines the dual oscillation mode with 
four parameters: amplitudes Ag and azimuthal wave 
numbers mr, mg. 


3. Investigated model of oscillating slender torus 

3.1. Location, thickness and frequencies identification 

Radial and vertical angular frequencies Dr -,^0 of the dual 
oscillation mode are given as a linear combination of eigen- 
frequencies ar,ag and Keplerian angular velocity LIko- 
Therefore the location of the torus centre in equilibrium rp, 
value of the spin of the central Kerr black hole and wave 
number pair mr, mg fully determine the ratio of oscilla¬ 
tion angular frequencies Qg/flr- Considering the properties 
of epicyclic and orbital frequencies in the Kerr spacetime 
(|T6rok fc Stuchlid l2005ll . our model identifies the lower 
and upper kHz QPOs frequencies with ui = |flr/27r| and 
Uu = \D.g/2'K\, respectively. As the aim of the article is to 
model twin peak HF QPOs with peaks frequency ratio close 
to 3 : 2, we choose rp in such a way, that the ratio Dig/Dr 
is just equal to 3 : 2 for a given wave number pair, In 
each of such positions of the torus centre we set the radial 
extent of the torus cross-section to rp/IO. Corresponding 
values of the parameter (3 are in accordance with slender 
torus condition 0. Arbitrary amplitudes of oscillations are 
fixed by setting Ar^ = Cdr,g- 

3.2. Optical properties of the torus 

The torus described above is assumed to be optically thick. 
Moreover, we use two other very simple assumptions. The 
torus surface emits radiation isotropically in its comoving 
local frame and the local flux integrated over the surface 
area of a thin vertical slice of the torus is conserved. Such 
a surface area is proportional to rc{tem,4>) ^ C{tem,4>), 
where C{tem,(j)) is the torus cross-section circumference, 
tern is the time of emission and rc{tem,</') is the radial 
coordinate of the centre of the torus cross-section. In the 
case of the used approximation, the investigated dual os¬ 
cillation mode describes pure radial and vertical displace¬ 
ments, and the torus cross-section circumference remains 
constant. Therefore, using normalisation to I for the equi¬ 
librium state, the local emitted intensity simply reads 

l-em {fern , (j)) = ro/rc{t em 1 4>). (33) 


® Investigated torus model does not consider the presence of 
the resonant coupling of oscillation modes, which can cause am¬ 
plification or excitation of QPOs on preferred values of radial 
coordinate (lHorakll2008l l. 
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3.3. Epicydic resonance axisymmetric setup 


The first investigated torus setup combines pure epicydic 
axisymmetric vertical and radial oscillation, where ra¬ 
dial and vertical oscillation frequencies are identical with 
radial i>r = Qr^Ko/'^'^ a-nd vertical ug = ujgElKo/‘2n 
epicydic frequencies, respectively. This setup corresponds 
to the often quoted epicydic resonan ce (ER) HF QPOs 
model (|Abramowicz fc KluznialJl200ltl based on the pres¬ 
ence of non-linear resona nt phenomena between ep i cydic 
disc oscillation modes (iKluzniak fc Abramowic i 120011 : 


_ _ s (IF ,, , - - - 

lAbramowicz et al.l l2003al lbl: iHorak 2008 1. The dual oscil¬ 
lation mode with such behaviour is related to wave number 
pair nir = 0, mg = 0. Therefore the frequency relation de¬ 
termining the radial coordinate of the torus centre rp reads 




me=0 


TTlr — O 


Ur 


(34) 


3.4. Relativistic precession-like non-axisymmetric setup 


The dual oscillation mode related to wave number pair 
mr = —1, mg = —2 yields the frequency relation in the 
form 

me = -2 


^ _ 2uk - I'e _ 3 

Vi vk -Vr 2 ’ 

mr = —1 


(35) 


where = Hxo/27r is the Keplerian orbital frequency 
at tq. As the denominator matches the periastron preces¬ 
sion frequency, in the Schwarzschild case the relation corre¬ 
sponds exactly to the relation of the relativistic precession 
(RP) QPOs model 


VI 


RP 


VK 


VK - Vr 


(36) 


The RP mode l was proposed in a series of papers by 
IStella fc Vietril (jl998lfl9^ 120021 1 i lMorsink fc Stelld (jlOOOl l 
and explains the QPOs as a direct manife station of rela¬ 
tivist ic epicydic motion of radiating blobs (IStella fc Vietril 
[19^ . In the case of slow rotation the frequency rela¬ 
tion (I35II still almost coincides with the ratio of twin 
peaks HP QPOs fr equencies predicted by the relation (l36)l 
(jTorbk et al.ll2012l l. Such an approach can be understood as 
a redefinition of the modulation mechanism of RP model, 
however, it preserves predictive power of the model. 


4. Numerical modelling of radiation emission, 
propagation and detection 

4.1. Ray-tracing in the Kerr spacetime 


Relativistic ray-tracing is a key ingredient of proper models 
of relativistic imaging, lightcurves of accretion structures, 
related power spectra as well as relativistic spectral line 
profiles. Different ray-tracing techniques were developed 
using direct numerical integra tion, transfer functions 
and el liptic integra l s (see , e.g.. ICunningham fc Bardeen 


Karas et al. 


Rauch fc Blandfor 



Viergutj~ 19931: Schnittman 


Beckwith fc Pond l2005t 
Broderick fc Loebll2005l : iBakala et ^120071: [Pexter fc Agoll 


20091: Vincent et al. _ 20111: Schee fc StuchlikI 120131 : 

Chan et al.1 l2013[ lYang fc Wangl l2014ll . Our parallel 


code LSDplus uses reverse ray-tracing implemented by a 
direct numerical integration of the null geodesics. 

Components of the four-momentum of a photon in the 
Kerr spacetime are given by 

= r = SrE“^Y^i?A.q(»'), 

= 0 = SgTi ^ 0A,q(0) , (’37^ 

p'i^ = (j) = [2ar -|- A — 2r) cosec^fl] , 

p* = i = (E^ _ 2arA) , 

where dotted quantities denote differentiation with respect 
to some affine parameter, and the sign pair Sr,sg describes 
orientatio n of radial and latitudinal evolut i on, respectively 
(see, e.g., ICarteil[l968l : iMisner et al.lll97^ IChandrasekharl 
Il983ll . Radial and latitudial effective potentials read 


R\,q (r) = [(r^ + a^) - aA]^ 


A 


q- {X 



©A,9 {d) = q-\- cos^ 9 — A^cot^0. 


(38) 


Here A, q are constants of motion related to the photons 
covariant angular and linear momenta. The LSDplus code 
performs a time-reverse integration of the package of null 
geodesics falling on a virtual detector of a distant observer 
with the screen resolution 1000 x 1000 pixels located at 
{9obs, r = 1000 M, (^ = 0) and traces the intersection of 
the geodesics with the disk surface corresponding to the 
emission events. 

The code LSDCode enables modelling relativistic opti¬ 
cal effects in the sky of the observer located anywhere above 
the event horizon of a Kerr black hole. In the local reference 
frame related to the such general observer, the components 
of the four-momentum of photons with energy normalized 
to one, falling on the pixel with coordinates a:, y, can be 
written as follows (see the right side of Figure [T]): 


kt = -1, kr = \J\- kg = -y, k^ = -x. (39) 


Then one can obtain the coordinate covariant components 
of the four-momentum by transforming the local compo¬ 
nents dSi), using appropriate frame tetrads of one-form by 
relation 

The frame tetrad of one-form related to static observer in 
the Kerr spacetime is given as 


AO _ 


^ 2r ^ ^ 2ar sin^ 6 


eft) = 

e(«) = 


{o,v/^,o,o}, 

{0,0,yE,0}. 


— 


0,0,0, A / ^ 


(41) 


The constants of motion X{x, y) and q{x, y) can be easily 
obtained by straightforward calculation from the compo- 
nent s of the four-moment um (I40L using the relations (see, 
e.g.. IChandrasekhailll983[ l 

A = (42) 

Pt 

+ (Atan(|-0))'-a2cos2 0. 
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Figure 2. Examples of the frequency shift maps for different distant views of the oscillating slender torus. The colour 
boxes on the right display the false colours scale of frequency shift values. The left column corresponds to the case of 
the ER axisymmetric dual mode. The right column corresponds to the case of the RP-like non-axisymmetric dual mode. 
The rows of the figure correspond to spin values a = 0, a = 0.5 and a = 0.96, respectively. 
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However, in the investigated case of a distant static observer 
{fobs —t oo), when the rays reaching the observer position 
are almost parallel (see Figure [T]), the relations between 
coordinates on the detector scre en and constants of mo¬ 
tion c an be simply written as le.g. lCunningham fc BardeenI 
1197.111 

X = —— , y = Q\,q (Oobs) ■ (43) 

sm Oobs 

In the event of detection, the constants X{x,y) and 
q{x,y), together with the initial conditions (coordinates 
of observer and sign pair Sr-,se) fully determinate the re¬ 
verse temporal evolution of the zero geodesics of photons 
falling on a pixel of the detector screen with the coordi¬ 
nates x,y. The used Runge-Ku tta method of the eighth or¬ 
der (Dorman-Prince method) ([Press et alJl200^ integrat¬ 
ing the null geodesics reaches the relative accuracy of 10“^®, 
which, in the case of a central black hole with stellar mass 
M = 5 Mq, corresponds to the order of accuracy 10“^^ me¬ 
ters on a radial coordinate. In order to determine the proper 
orientation of the radial and latitudinal component of the 
four-momentum (EZl), the code additionally analyses the 
positions of the radial and vertical turning point and sets 
the corresponding signs Sr and sg. The integration of equa¬ 
tions (1571) by the the Runge-Kutta method of the eighth 
order proceeds naturally with adaptive step. However, the 
resulting null geodesic is then finally interpolated by the 
polynomials of the 3rd order for the chosen equidistant time 
step AT, in the case of central black hole mass M = 5 Mq, 
corresponding to 10“^ sec . 


the event of emission on the surface of the torus. lem is local 
intensity on the surface in the comoving frame given by the 
equation (l33l) . The total instantaneous detected bolometric 
flux F{tobs) is a sum of intensities (1451) detected by individ¬ 
ual pixels multiplied by the solid angle AH subtended by 
the pixel in the observer sky 

1000 1000 

F{tobs) = EE /y (tobs) ah, (46) 

i=i j=i 

where i is the index of a pixel column and where j is the 
index of a pixel ro’wQ- Since we are using relative units, 
we can simply set AH = I. The resulting lightcurves are 
generated in the time resolution AT corresponding to 20 
time samples per characteristic vertical oscillation period 
of the analyzed dual oscillation mode. 

4.3. Power spectra 

In order to calculate power spectral densities (PSD), the 
resulting lightcurves are processed by fast Fourier trans¬ 
form (FFT). The power spectral density at frequency fk = 
k/{NAt) is given as the square of the modulus of the FFT 
of the signal as 


PSDifk) 


1 

N 


N-l 

E -P'(^i)exp(27rijA:/A) 


i=o 


2 


(47) 


4.2. Radiating surfaces and lightcurves 

The surfaces of the tori are modelled by the grid with res¬ 
olution 15 nodes in the radial direction and 75 nodes in 
the azimuthal one. The grid contains the time-dependent 
information about coordinates, local intensity, and the four- 
velocity in the nodes. The time resolution of the surface of 
the tori necessarily corresponds to the time resolution of 
the interpolation steps of the geodesics package AT. The 
code LSDplus traces the intersections of the geodesics and 
linearly interpolated surface of a torus between triads of the 
nodes of the grid. Assuming the normalized energy in the 
local observer’s frame, the frequency ratio of the emitted 
and observed radiation from the torus can be expressed us¬ 
ing projection of the four-momentum of a photon to the 
four-velocity of the surface of the torus in the event of 
the emission as follows: 


where F{tj) is the observed flux (1461) at observation time 
tj = j At and N is the total number of time samples in the 
lightcurve. 

4.4. Iron Ka line profiles 

The fluorescent iron Ka line consists of two comp onents 
with FWHM Ri 3.5 eP and separation r; 13 eP fsee iBask^ 
I1978L for details). Considering the energy resolution AE = 
10 eV of the simulation code, we approximated the rest the 
rest iron Ka line profile by Lorentzian peak with the scale 
factor 7 = 20 eV located on Eq = 6.4fcey. The instanta¬ 
neous observed flux per pixel in the energy bin with the 
central energy Ec is given as 

^{tobs,Ea.) = lemitobs “ tdelay)g^f {Ec/g, 7, Eq) AH, (48) 

where Lorentzian function / reads 


Radial and vertical components of the surface four-velocity 
are given by equations (I25I29I) . The remaining components 
u™’', can be easily obtained using normalization con¬ 
dition = — 1 together with the assumption of constant 
specific angular momentum (IVincent et al.l[20T3l . Then the 
instantaneous bolometric intensity detected by each pixel 
of the screen of a small virtual detector is calculated as 

fobsifobs) — f-em {fobs l'delay)g j (4b) 

where t^eiay corresponds to the time delay (the change of 
a time coordinate) along the appropriate photon trajectory 
connecting the event of detection of a photon on a pixel with 


f{x,^,xo) 


1 

^7[l + (“)^]’ 


(49) 


Then the observed instantaneous iron Ka line profile is 
constructed by summing energy bin fluxes per pixels (1481) 
over all pixels of the virtual detector in the given time sam¬ 
ple. The final integrated line profile is obtained by summing 
instantaneous line profiles over all time samples. 


^ We assume tiny angular size of the torus image in the ob¬ 
server sky and therefore constant AH. 
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5. Methods 

Applying ray-tracing methods described in the previous 
section, we performed simulations of lightcurves of the 
above discussed model of a slender accretion torus oscil¬ 
lating in the dual mode regime. In order to obtain closed 
trajectories of the torus surface and corresponding closed 
lightcurves, we simulated behaviour of emitted radiation 
during three periods of vertical oscillation considering the 
hxed ratio of oscillation frequencies Vujvi = 3/2. Then we 
calculated power spectra of obtained lightcurves by rela¬ 
tion dUl). We studied the impact of the central black hole 
spin using its three representative values a G (0, 0.5, 0.96). 
The upper limit of the investigated spin value is almost 
the highest one, for which the location of equilibrium torus 
centre vq can be found by the RP-like frequency relation 
((551) . The impact of the observer inclination (polar angle) i 
is analyzed using the following set of inclination values: 

TT TT TT TT 5 5.5 5.75 TT \ 

, •, —1 —1 —1 TT, TT, -TT, — 1 

’ 12’ 6 4’ 3 12 ’ 12 ’ 12 ’2/ 

We analyze the magnitude relations of five prominent PSD 
peaks corresponding to two main dual mode oscillation fre¬ 
quencies, their higher harmonics and their sum or differ¬ 
ence, particularly summarized in the Table [TJ Power spec- 


Frequency scaled in i/j 

Peak origin 

1/2 

I'u — Vl 

1 

VI = |Dr/27r 

3/2 

Vu = Ds/27r 

2 

2vi 

5/2 

Vu -f Vl 

3 

, 2i^u, 


Table 1. The most prominent peaks observed in power 
density spectra of simulated lightcurves. 


tra are calculated for all combination of investigated values 
of the spin and the inclination. We plot magnitudes of these 
peaks as linearly interpolated functions of the inclination 
for all investigated values of the spin and for both investi¬ 
gated torus setups separately (see left panels of Figures [3] 
and ID). 

Moreover, the integrated iron Ka line profiles emit¬ 
ted from torus surface are modelled for inclination values 
* ^ f’ f’ f’ §)■ geometry of the torus optical 
projection is illustrated in Figure [5] by the frequency shift 
maps on the virtual detector screen drawn in false colours 
related to the frequency shift values 0. The maps clearly 
show that the simulation parameters are sufficient to distin¬ 
guish the first three relativistic images. Also computed iron 
Ka line profiles displays significant secondary blueshifted 
horns related to secondary (first indirect) relativistic im¬ 
ages (see right panels of Figures 15151) . 

6. Results 

6.1. Results for ER axisymmetric setup 

Table [5] summarizes parameters of the slender torus oscillat¬ 
ing in the pure epicyclic axisymmetric dual mode (m^ = 0, 

® The frequency shift maps displayed in Figure [2] are not iden¬ 
tically zoomed and positioned with respect to the whole observer 
sky. 


Spin 

0.0 

0.5 

0.96 

P 

0.04 

0.04 

0.04 

Central radial coordinate ro 

10.80 

7.92 

4.35 

Torus radial extent 

1.08 

0.79 

0.44 

Torus vertical extent 

0.80 

0.61 

0.38 

Max. radial displacement 

0.27 

0.20 

0.11 

Max. vertical displacement 

0.20 

0.15 

0.09 


Table 2. Values of parameters describing the investigated 
slender torus model for ER axisymmetric setup. 


m 0 = 0) for three investigated values of the spin. The radial 
coordinate of equilibrium torus centre tq given by relation 
(1551) remains above the black hole photosphere^ as well as 
the ergosphere for all such configurations. 

The lightcurve waveforms are influenced by a complex 
interplay of the general relativistic frequency shift, the time 
variation of the emitting torus surface area and the time 
variatio n of the apparent to rus area on the virtual detector 
screen ((Bakala et al.l [201511 . All these effects strongly de¬ 
pend on the radial coordinate of the emission event, the in¬ 
clination of the distant observer i and the central black hole 
spin a. Nevertheless, the magnitude relations depicted on 
left panels of Figure [3] exhibit certain identical qualitative 
features for all investigated spin values. The pair of peaks 
corresponding to radial vi and vertical torus oscillation 
frequency is dominant on wide range of i. For this inclina¬ 
tion range the examined ER axisymmetric setup predicts 
twin peaks HE QPOs frequency ratio equal to 3/2 in ac¬ 
cordance with expectations. The ui peak remains the most 
prominent for small and medium values of distant observer 
inclination i. The magnitude of the peak grows with i 
becoming the most prominent for relatively high inclina¬ 
tions, but it rapidly falls down for exact or almost exact 
equatorial observers. At the same time, the magnitude of 
the peak corresponding to 3ui, 2vu rapidly grows. Therefore, 
in the case of such observers, the examined torus configu¬ 
ration predicts frequency ratio equal to 1/3 for the pair 
of the most distinguishable HE QPO peaks. In the case of 
zero or moderate spin (a = 0, a = 0.5) and small inclination 
{i < ^)) the magnitude of the 2 ui peak slightly exceeds the 
ui peak magnitude and the predicted twin peaks HE QPOs 
frequency ratio value is equal to 1/2. As illustrated on the 
plots of magnitude relations (see left panels of Eigure[3]), 
the spectral content of less distinct PSD peaks varies de¬ 
pending on values of i and a. 

Like the power spectra behaviour, the Iron Ka line pro¬ 
files keep some identical qualitative features for all investi¬ 
gated spin values. Naturally, the energy span of line profiles 
is expanded and shifted down with decreasing tq. The de¬ 
pendence of the height of the primary blueshifted horns on 
i is the only qualitative difference observable for different 
spin values (see right panels of Figure [3]). 

6.2. Results for RP-like non-axisymmetric setup 

Table [3] summarizes parameters of the slender torus oscillat¬ 
ing in the RP-like non-axisymmetric dual mode {rrir = — 1, 


° Kerr black hole photosphere - a region of spherical unstable 
photon orbits reaches maximum extent in the equatorial plane 
between corotating and counter-rotating circular photon orbits, 
while it b ecomes infinitesimally thin on the polar axis (see, e.g., 
[T^[2003. for details) 


* G 0.01 
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Figures. The simulations outputs in the case of the ER axisymmetric dual mode. Left panels: The amplitudes of 
prominent PSD peaks (see Table [1]) as a function of distant observer inclination i . Right panels: Iron Ka line profiles 
constructed forfs The rows of the picture correspond to the cases of spin values a = 0,a = 0.5 and 

a = 0.96, respectively (see Tabled). 


mg = —2) for three investigated values of the spin. In the 
case of zero or moderate spin (a = 0, a = 0.5), the ra¬ 
dial coordinate of equilibrium torus centre rg given by rela¬ 
tion (1351) remains located above the black hole photosphere 
as well as the ergosphere, while in the case of high spin 
(a = 0.96), the torus is located inside the ergopshere and 
therefore also deeply inside the photosphere. Moreover, in 


the case of high spin, the radial extent of the torus obtained 
as rg/lO exceeds the location of the cusp of equipotential 
surfaces (|Blaes et al.ll2006HStraub &: Sramkoval2009ll . Such 
a torus configuration becomes unstable. Unfortunately, a 
surface area of the high spin stable configuration with max¬ 
imum possible 13 = 7.0.10“^ is almost negligible and emit¬ 
ted flux is comparable to numerical error of the simula- 
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Figure 4. The simulations outputs in the case of the RP-like non-axisymmetric dual mode. Left panels: The amplitudes 
of prominent PSD peaks (see Table [T|) as a function of distant observer inclination i . Right panels: Iron Ka line profiles 
constructed foriS ^). The rows of the picture correspond to the cases of spin values a = 0,a = 0.5 and 

a = 0.96, respectively (see Table [3|). 


tion. Therefore, we keep the torus radial extent equal to 
ro/10 and choose the unstable high spin configuration with 
P = 0 . 002 . 

In the case of zero or moderate spin (a = 0, a = 0.5), the 
qualitative picture of power spectra and iron Ka line pro¬ 
files behaviour is very similar to the case of the ER axisym- 
metric setup discussed in the previous section (see Figure 
a. The prediction for low inclination represents the main 


difference. In the Schwarzschild case, the pair of peaks cor¬ 
responding to radial i^i and vertical oscillation frequency 
remains dominant also for * < ^ predicting the twin peaks 
HF QPOs frequency ratio equal to 3/2, as depicted on the 
top left panel of Figure 0] In the case of a = 0.5 and a very 
small inclination, the magnitude of the — i^i peak (in¬ 
stead of the 2iyi peak acting the same way in the case of the 
ER axisymmetric setup) slightly exceeds the z/; peak mag- 
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Spin 

0.0 

0.5 

0.96 

P 

0.02 

0.02 

0.002 

Central radial coordinate ro 

6.75 

4.66 

1.85 

Torus radial extent 

0.68 

0.47 

0.18 

Torus vertical extent 

0.27 

0.18 

0.02 

Max. radial displacement 

0.17 

0.12 

0.05 

Max. vertical displacement 

0.07 

0.05 

0.004 


Table 3. Values of parameters describing the investigated 
slender torus model for RP-like non-axisymmetric setup. 


nitude, and the predicted twin peaks HF QPOs frequency 
ratio is equal to 1/2 (see the middle left panel of Figure H]). 

The unstable high spin configuration exhibits an en¬ 
tirely different picture, as shown in bottom panels of Figure 
01 In the whole range of inclination, the pair of dominant 
PSD peaks corresponds to the radial vi frequency and its 
second harmonic 2r'/. Therefore, in the unstable high spin 
case, the examined RP-like non-axisymmetric setup sur¬ 
prisingly predicts the twin peaks HF QPOs frequency ratio 
value equal to 1/2 only. The dramatic change of the opti¬ 
cal projection properties is documented by iron Ka line 
profiles in the bottom right panel of Figure U) The en¬ 
ergy span is significantly shifted down, and both primary 
and secondary blueshifted horns have comparable heights. 
The computed profiles also display numerical instabilities 
caused by small emitting surface area of the unstable high 
spin torus configuration. The line profile related to i = ^ 
exhibits exceptional behaviour, as its extremely wide en¬ 
ergy span reaches 11 keV and its primary blueshifted horn 
is absolutely dominant. Corresponding frequency shift map 
in the bottom right panel of Figure [2] also illustrates the 
equatorial optical projection of the torus located in the 
close vicinity of the Kerr black hole event horizon. It is 
clearly visible, that the angular size of both secondary and 
tertiary relativistic images exceeds the angular size of the 
primary image. Moreover, despite the signihcant gravita¬ 
tional redshift, the high Keplerian orbital velocity causes 
high Doppler blueshift of the left side of the torus projec¬ 
tion corresponding to the extreme height of the primary 
blueshifted horn in the iron Ka line profile. 

7. Conclusions and perspectives 

The aim of this article is to model twin peaks HF QPOs 
as a spectral impact of isotropically radiating slender tori 
oscillating in a dual mode regime i n the close vic i nity o f 
Kerr black holes. It was shown by I Vincent et al.l (|2014ll 
that signihcant fraction of the observed flux is regulated 
by the torus motion described by the lowest oscillation 
modes. Therefore, we examined two conhgurations of the 
dual oscillation regime based on the lowest radial and 
vertical oscillation modes with different azimuthal wave 
numbers. Appropriate frequency relations correspond to 
the two competing QP Os models, the epicyclic reson ance 
(ER) HF QPOs model (lAbramowicz fc Klulniakll2001f) and 
slightly modihed rela tivistic precession (RP) QPOs model 
(|Stella fc Vietril[r999ll . We model twin peaks HF QPOs by 
the pair of the most prominent peaks in the obtained power 
spectra. Our results show that independently of the spin, 
the ER axisymmetric setup yields power spectra with the 
pair of dominant PSD peaks corresponding to the frequen¬ 
cies of radial and vertical oscillation modes with proper 
ratio equal to 3/2, except some special cases of very high 


or very low distant observer inclinations, where higher har¬ 
monics becomes prominent. The predictions of the RP-like 
non-axisymmetric setup are almost identical to the ER case 
for zero or moderate spin configurations. An unstable high 
spin RP-like configuration with the torus located in the er- 
gosphere exhibits dominant PSD peaks pair corresponding 
to the frequency of radial oscillation modes and its second 
harmonics in the whole range of inclinations and therefore 
it predicts constant frequency ratio equal to 1/2. The entire 
change of the optical projection is also documented by the 
different iron Ka line profiles with respect to the previous 
cases. 

The analysis presented in the article primarily focuses 
on relative ratios of the amplitudes of the most prominent 
frequency peaks in the modelled power spectra. However, 
the absolute fractional rms of the amplitudes of HF QPOs 
peaks in the detected signal depends not only on the am¬ 
plitudes of perturbation A^^g, but also on relations of in- 
dividual componen ts of the whole global source flux. In 
iBakala et al.l (|2014ll we defined an empirical model of global 
source flux, which mimics the so-called high steep power 
law (HSPL) state in GRS 1905-1-105, including steep spec¬ 
trum and power-law dominated variability with an ad- 
ditional broad Lorentzian com ponent at low frequencies 
(jMcClintock fc Remillardll2006f) . Then we used this back¬ 
ground for the analysis of resolution of HF QPO peaks 
for a similar but simpler model of oscillating slend er torus 
in the Schwarzschild geometry (iBursa et al.ll2003) . slowly 
passing the resonant orbit vq. Considering the capabilities 
of the RXTE and LOFT instruments simulated by their 
response matrices, we have shown that the present avail¬ 
able observational technology enables good detection of 
the pair of the most prominent peaks in such a modelled 
signal. Nevertheless, the behaviour of spectral content of 
less distinct PSD peaks remains probably the key task for 
the data analysis of future sensitive space observatory mis- 
sions for X-ray timing, e.g. the proposed LOFT mission (see 
iBakala et al.H20l4lFeroci et al.ll201^lKaras et al.ll20il. for 

details). 

Our simulation yields iron Ka line profiles which are 
very different from the line profile integrated over the en¬ 
tire accretion disk. Therefore, the next related topic for 
future space X-ray missions can be sensitive frequency- 
resolved spectroscopy, which will be able to isolate the 
spec tral component oscillating on the QPO freque ncy (see, 
e.g.. lAxelsson et ^l20l4 iRevnivtsev et al.lll999li . Phase- 
resolved spectroscopy tracing the iron line profile changes 
throughout an oscillation cycle could be another poten¬ 
tially interesting diagnostic tool. Such idea has already 
been suggested for the case of low frequency QPOs (e.g., 
lingram fc Donel[2012bl: [Tsang fc Butskvll2M^ . but its ap¬ 
plication for the studied model of HF QPO driven by slen¬ 
der torus oscillations will require a more detailed future 
study. 

We assume an optically thick slender torus, and thus 
the influence of the torus opacity can be also subject for 
our future research. It is possible to generalize the exam¬ 
ined slender torus model considering p ressure effects on the 
mode frequencies and the torus shape (IStraub fc Sramkoval 
I 2 OO 9 II . In the presented study, the location of the torus cen¬ 
tre is empirically chosen using the observed twin peaks HF 
QPOs ratio. The resonant coupling of oscillation modes can 
prefer particular values of ra dial coordina te for amplifica¬ 
tion or excitation of QPOs (jHorakI [2008fl . Therefore, fur- 
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ther important improvement of the simulations can incor¬ 
porate mode resonant coupling. Future implementation of 
improved models of non-slender tori into the used LSDplus 
code can be the next step towards more realistic results. 
A future detailed comparative study of spectral harmonic 
content can also be devoted to the possibility to distinguish 
between QPOs models based on either accretion tori oscil¬ 
lations or orbital motion of radiating blobs. 
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